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Summary. - A classical model of the hydrogen atom in a static electric 
field is studied, basing upon the work [ Hooker A. et al, Phys. Rev. A, 55 
(1997) 4609 ]. In that work the electrons are supposed to move along Kepler 
orbits around the nucleus, while interacting with the external field. That 
classical model reproduces very well the true energy shift of the linear Stark 
effect. The agreement with the second order effect is poor. It is shown here 
that the results for the quadratic Stark effect may be considerably improved 
if the electrons are still allowed to move along classical trajectories, but whose 
initial conditions are statistically sampled from a distribution resembling the 
quantum mechanical one. 
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PACS 03.20.+i - Classical mechanics of discrete systems: general mathemat- 
ical aspects 
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1. - Introduction 



The use of classical models to study quantum mechanical systems is largely 
employed in microscopic physics. For large quantum numbers it is justified on 
the basis of the correspondence principle. Also at small quantum numbers its 
utility is not negligible since it allows at least to obtain a qualitative insight 
of the dynamics of the system; in the best cases, it is even possible to extract 
quantitative results in a easy way and with a good accuracy. Some examples 
are the scattering of electrons from an atom [[I], |2[], and the treatment of 
energetic ion-atom collisions 0, |], |5| . Classical model are obviously favoured 
over their quantum-mechanical counterparts from the numerical point of 
view. In those cases where also the quantum problem is solvable, the classical 
picture is important since it allows to clearly visualize the problem at hand. 
For these reasons the correspondence between the classical and quantum 
mechanical descriptions is currently a well studied topic. 

Quite recently it has been proposed a purely classical view to the Stark 
effect in hydrogen ffl : the field is supposed to interact with classical electrons 
following Kepler orbits around the nucleus. Following this work, in ref. the 
energy shifts for the linear and quadratic Stark effect have been computed. 
The agreement between the classical result and the exact quantum mechan- 
ical one is excellent for the linear effect but is not good for the quadratic 
effect, which only asymptotically for large quantum numbers approaches the 
correct value. 

The purpose of the present paper is to show that, within a purely classi- 
cal formalismf] very similar to that of [0, the quadratic Stark effect may be 
reproduced very accurately: it is enough to relax the bounds on the electron 
trajectory and adopt a statistical approach: here, the electron is still mod- 
elled as a classical particle following Kepler orbits but the initial conditions 

2 It must be stressed that any " classical" model necessarily incorporates some elements 
of the quantum theory in order to simulate a microscopic object. The relevance of the 
approach Q and of the present paper is that only initial conditions are to be chosen 
compatible with the laws of quantum mechanics. 
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are picked up from a distribution subject to certain rules. It will be shown 
that a much better agreement is obtained by this slightly more sophisticated 
approach. 

Atomic units will be used throughout this work unless explicitly stated. 

2. - The classical model of [j7j] 

As well known, the study of the Stark effect is easily performed in parabolic 
coordinates, where the hydrogen atom is classified by the quantum num- 
bers n,rii,n2,m. We are considering small fields, for which a perturbative 
approach is adequate and they still are good quantum numbers. 

I will first briefly summarize the treatment of the first-order Stark effect, 
which will be useful also to introduce the main concepts used later. First of 
all, let us recall some results stated in |J. In that work a correspondence 
has been done between the quantum mechanical operators and their classical 
counterparts: 

• to each set of values n 1; n, m there corresponds a set of classical Ke- 
plerian orbits of the electron which are generally elliptical, with eccen- 
tricity s = a/1 — l 2 /n 2 (I is angular quantum number in spherical co- 
ordinates, and the classical angular momentum). When time-averaged 
over a period of rotation, the ellipse yields a non zero mean electric 
dipole moment: 

i r T 

< d >= - / zdt = 3/2n 2 e . (1) 
T Jo 

• The Lenz vector A = pxL — r/|r| is classically a conserved quantity, 
as is its ^-component in quantum mechanics: 

K = (ni - n 2 )n . (2) 

A is related to the electron orbit by |A| — e. Since, besides this, in 
it was shown that d and A point towards the same direction, one may 
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identify them through 

d = (3/2)n 2 A . (3) 

• The energy shift in presence of an electric field F is calculated from 

A£ (1) = -d F . (4) 

It is immediate to see that this result agrees with the energy shift 
obtained using quantum mechanics provided we use for d its form 
(D with A given by (|): 

AE^ = -^n( ni -n 2 )\F\ , (5) 

In it was verified by numerically integrating the Kepler orbits in pres- 
ence of an electric field and with initial conditions compatible with the quan- 
tized value (0) that the classical energy shifts from its unperturbed value of 
the quantity @. 

The quadratic effect is of relevance when rii = n 2 , in which case the 
linear term vanishes. This means-in the language of 0-that the mean dipole 
moment vanishes, and this happens when the orbits are circular and lying 
on a plane perpendicular to the field axis. The electric field may induce a 
dipole moment by shifting the electron and the nucleus with respect to each 
other. The energy shift for the quadratic Stark effect is defined by 

AE^ = -^aF 2 , (6) 

a being the polarizability. a is related to the induced dipole moment by 

aF = d . (7) 

The quantum mechanical value of a for states with n\ = n 2 and m = n — 1 
is@ 

4 

71 

O-qm = ~~r (4ri 2 + 9n + 5) . (8) 
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In [[7] a simple approximation is used to compute the classical value of a, 
a c f. be r the radius of the circular orbit of the electron, and 8z the shift 
along the direction of F between the electron and the nucleus induced by the 
external field. For small values of 5z we may approximate the Coulomb force 
on the electron as 5z/r 3 , which balances the force exerted on the electron by 
F when 

F = -z ■ W 

The induced average dipole moment is d = 5z and, from eq. (|SJ) and the 
definition of a (|7|), 

a d = r 3 = n 6 (10) 

since for circular Bohr orbits r = n 2 . 

Even if this is a very simplified model, it was verified in @ that it de- 
scribes very accurately the classical system: the Hamilton's equations for the 
electron were numerically solved in presence of the electric field. The electron 
energy was determined from its position and momentum. The numerical re- 
sults were found to agree well with eq. fllOD (see fig. 3 of 0). It is therefore 
correct to assume that ot c i is well approximated by eq. fllOf) . 

Eqns. (D and flj3) only agree for n — > oo; the greater discrepancies are 
for small ns: for example, when n — 1 the ratio is a c i/aQM = 2/9 ~ 0.22. 

3. - Improvements over the simple model 

Both quantum and classical mechanics admit a remarkable unified description 
in terms of fluid dynamics: the Schrodinger equation 

^^=(-^V 2 + ^(r,t))^(r,t) (11) 

through the replacement 

i[>(r,t) = ^p(r,t)exp Us^tfj (12) 
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may be rewritten into the set of two equations 

dp 



m + V • (pv) = (13) 
= -W e// (14) 



where 



U eff = U(r,t)-^^f = U(r,t) + U q (r,t) (15) 

v(r,t) = (yS)/m (16) 

The time evolution of a classical phase space distribution is given by the 
Liouville equation 

I = w < 17 > 

Lit) = V r U-V p - — -V r (18) 

m 

from which, multiplying by 1 and p respectively and then integrating, 

f + V-( P v)=0 (19) 

= -V r U - -V r Ii (20) 
dt p 

p(r,t) = J dpf(r,p,t) (21) 

v = — / dpp/(r,p,t) (22) 

U i:i = — dppjPjf(r, p, i) - pviVj (23) 

771 J 

The only difference between the two sets of equations is that the role of the 
quantum potential U q has been replaced by the stress tensor II. 

From / one gets the projection over the position and momentum coordi- 
nates 

p(r,t) = J dpf(r,p,t) (24) 
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p(p,t) = I drf(r,p,t) (25) 

In classical approximations the quantum dynamics of a system is com- 
puted through an averaging over an ensemble of classical trajectories. The 
average is done over the statistical distribution of initial conditions /, with 
/ chosen such as to closely reproduce the wave function in position or mo- 
mentum space. 

If / is to be stationary (as is the case here), it can depend only on con- 
stants of the motion ||. A largely employed choice is to make / depending 
only upon the energy: /(r, p, £) = f(E). There is not an unique choice pos- 
sible for / since-classically-does not exist any / such that both p and p are 
equal to their quantum mechanical counterpart. The choice is done on the 
basis of mathematical simplicity and accuracy. 

Once a choice for /-and therefore for p, p-is done, the calculations of the 
previous section still hold: simply, they are to be repeated for an ensemble of 
electrons, each of them is still obeying Newton's equations, and in particular 



follows a circular orbit. Therefore, eq. (|10|) is replaced by its average 



a=<r 3 > = J drr 3 p(v) . (26) 

Notice that the model of || [7j may be seen as a particular case, by putting 
p ~ S(r(t) —Tj((t)) , p ~ <Kp(£) ~ Pif(0)> whh r K , p K position and momentum 
corresponding to a Kepler orbit. 
I will consider three special cases: 

1. first of all, the true quantum mechanical distribution is used for p 

p(r) = |^ m (r)| 2 (27) 

with I = \m\ — n — 1. This is an obvious choice, and is done in order to 
give an insight of the effectiveness of the method when applied under 
those which should be the conditions closest to the true ones. It has 
the defect that p now encompasses a region classically forbidden to the 
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electron: from 

p 2 1 1 



(28) 



2 r 2n 2 
one gets r < 2n 2 . 

2. The second choice is therefore to use a "truncated" distribution: 

jc\^ nlm (r)\ 2 ,r<2n 2 

P = S , ( 29 ) 

10 ,r>2n 2 

where C stands for a normalization constant. 

3. Finally, I use a microcanonical distribution 

^sM B+ ^) (30) 

from which one obtains, after substitution in eq. 

i o 1/2 



"H;-s?J (31) 

which is a real quantity only for r < 2n 2 . This microcanonical distribu- 
tion is largely used in classical calculations of ion-hydrogen scattering 
processes for generating initial electron distributions, since it has the 
property of correctly reproducing the true electron momentum distri- 
bution p ||. 

All the integrals fl26|) may be performed analytically using the ps of items 
1-3. In table 1 I report the values for the polarizability obtained using 1-3, 
together with the values from eqns. (0) and ([TO]) . It is clearly discernible 



the slow convergence of a c i to cxqm if compared, in particular, to «2 which 
already at n = 1 is within 20% of the true value. It is curious, on the 
other hand, the behaviour of a 3 which is always very close to oqm and, for 
n = 4, 5 is the best approximation, but becomes worse at large us. Another 
point to remark is that, for large n, ax and ct2 become nearly equal, as it 
should be from their definition. On the average, using any of the suggested 
distributions allows at least to halve the error with respect to the results of 
ref. 0. 
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4. - Conclusions 



In ref. [0] it is remarked that the main motivation of their work was to 
show that some quantal results may be recovered, at least partially, within 
a purely classical framework, provided only that the initial conditions be 
chosen as compatible as possible with the laws of quantum mechanics. This 
paper follows exactly that line of thought, being a refinement of that work in 
that the initial conditions have been chosen in a more correct manner, but 
still remaining within a classical description of the system. One may wonder 
if the price paid to have this greater accuracy is too high, since we used 
exact wave-functions, which is equivalent to solve the quantum mechanical 
problem. It is not so, since significant improvements are obtained using no 
matter which distribution, even that of eq. ([H]) which is based upon purely 
classical considerations. The other two distributions have been chosen just 
to provide the reader with a comparison. This is important for all these cases 
where extracting the wave function is too difficult and one is forced to resort 
to approximations. As already pointed out in [0 , it could be worth exploring 
in some of such situations. 

This work has been supported by a grant of the Italian MURST. The hospi- 
tality offered by the Consorzio RFX is acknowledged. 
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n 




0L 2 


a 3 




OlQM 


1 


15/2 


3.616 


8/5 


1 


9/2 




(1.67) 


(0.80) 


(0.36) 


(0.22) 




2 


210 


127.6 


512/5 


64 


156 




(1.35) 


(0.82) 


(0.66) 


(0.41) 




3 


1701 


1350.6 


5832/5 


729 


1377 




(1.24) 


(0.98) 


(0.85) 


(0.53) 




4 


7920 


7069.7 


32768/5 


4096 


6720 




(1.18) 


(1.05) 


(0.98) 


(0.61) 




5 


53625/2 


25313 


25000 


15625 


46875/2 




(1.14) 


(1.08) 


(1.07) 


(0.67) 




6 


73710 


71561 


373248/5 


46656 


65772 




(1.12) 


(1.09) 


(1.13) 


(0.71) 




10 


1.328 10 6 


1.325 10 6 


1.6 10 6 


1 10 6 


1.2375 10 6 




(1.07) 


(1.07) 


(1.29) 


(0.81) 





Table 1: polarizability calculated from the different methods explained 
in the text. 01,2,3 from choices 1, 2, 3 respectively; a c i from eq. (0); cxqm 
quantum mechanical value from eq. (|]). Between parentheses are reported 
the ratios ajaqu- n is the principal quantum number. 
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